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ABSTRACT 



We discuss the effect of the finite size of galaxies on estimating small-scale 
relative pairwise peculiar velocity dispersions from the cosmic virial theorem 
(CVT). Specifically we evaluate the effect by incorporating the finite core radius 
Tc in the two-point correlation function of mass, i.e. ip{r) oc (r + Tc)'"' and the 
effective gravitational force softening on small scales. We analytically obtain 
the lowest-order correction term for 7 < 2 which is in quantitative agreement 
with the full numerical evaluation. With a nonzero and/or Tc the cosmic 
virial theorem is no longer limited to the case of 7 < 2. We present accurate 
fitting formulae for the CVT predicted pairwise velocity dispersion for the case 
of 7 > 2. Compared with the idealistic point-mass approximation (r^ = Tc = 0), 
the finite size effect can significantly reduce the small-scale velocity dispersions 
of galaxies at scales much larger than and r,,- Even without considering the 
finite size of galaxies, nonzero values for Tc are generally expected, for instance, 
for cold dark matter (CDM) models with a scale-invariant primordial spectrum. 
For these CDM models, a reasonable force softening Vg < 100/i~^kpc would 
have rather tiny effect. We present the CVT predictions for the small-scale 
pairwise velocity dispersion in the CDM models normalized by the COBE 
observation. The implication of our results for confrontation of observations of 
galaxy pair-wise velocity dispersions and theoretical predictions of the CVT is 
also discussed. 



Subject headings: cosmology: theory — large-scale structure of the universe — 
methods: statistical 
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1. INTRODUCTION 

The three-dimensional distribution of galaxies observed in redshift surveys differs from 
the true one since the distance to each galaxy cannot be determined by its redshift z only; 
for 2; <^ 1 the peculiar velocity of galaxies, typically ~ (100 — 1000)km/sec, contaminates 
the true recession velocity of the Hubble flow (e.g., Davis & Peebles 1983; Kaiser 1987; 
Hamilton 1992), while the true distance for objects at z ^ 1 sensitively depends on 
the (unknown and thus assumed) cosmological parameters. This hampers the effort to 
understand the true distribution of large-scale structure of the universe. Nevertheless 
such redshift-space distortion effects are quite useful since through the detailed theoretical 
modeling, one can derive the peculiar velocity dispersions of galaxies as a function of 
separation (Davis & Peebles 1983), and also can infer the cosmological density parameter 
Qo and the dimensionless cosmological constant Aq, for instance (Ballinger, Peacock & 
Heavens 1996; Matsubara & Suto 1996; Suto & Matsubara 1996). 

Theoretically, a conventional tool to predict small-scale relative pairwise peculiar 
velocity dispersions is the cosmic virial theorem (Peebles 1976; hereafter CVT). In its 
simplest form, the prediction is based on the idealistic assumption that galaxies are treated 
as point particles. As Peebles (1976) remarked, however, the finite size effect is very 
significant on scale ~ l/i~^Mpc even if one considers a radius of (20 ~ 100)/i~^kpc for 
typical galactic halos. In this paper we explore this realistic effect in more details and 
discuss its important impact on the application of the CVT to real observations. We will 
also give rehable predictions for the small-scale velocity dispersions of galaxies in COBE 
normalized CDM models. 

The plan of the paper is as follows; section 2 brieffy describes the CVT both in its 
simplest and in its more realistic form. The finite size effect is considered in §3, first based 
on a perturbation theory (§3.1) and then on a full numerical analysis (§3.2). In §4 we 
predict velocity dispersion of galaxies on small scales for COBE normalized CDM models. 
Discussion and implications of our result, and our main conclusions are given in §5. 
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2. COSMIC VIRIAL THEOREM 

Assuming that galaxy pairs are gravitationally in "statistical equilibrium" on small 
scales, Peebles (1976, 1980) derived an expression for their relative (one-dimensional) 
peculiar velocity dispersion as a function of their separation r: 

W) 

where p is the mean density of the universe, and ^p(r) and C,p are the two- and three-point 
correlation functions of mass. 



. 2 r w SGp f°° dr' f v' ■ 7. 



The observed two- and three-point correlation functions of galaxies, and ^g, are well 
approximated by the following forms (Groth & Peebles 1977; Davis & Peebles 1983): 

w = (2) 



r 

Cg(ri, r2, ra) = Qg [Uri)Ur2) + Ur2)Ur3) + UrsMri)] (3) 

where Tq = (5.4 ± 0.3)/i-^Mpc, 7 = 1.77 ± 0.04, Qg = 1.29 ± 0.21 {h is the dimensionless 
Hubble constant Hq in units of 100 kms~^ Mpc~^). Thus it is reasonable to assume that 
the two- and three-point correlation functions of mass, and (p, also obey the same scaling 
except for the overall amplitudes: 

= 4^gW' Cp(n, r2, ra) = -^Cgin, r2, rg). (4) 

Then a straightforward computation yields 

{vl2ir))cvT = 6GpQpr2r'-'^Io{l)/bl, (5) 

where 

^(7) 



(7-l)(2-7)(4-7) 
°° 1 + X " 



X / dx , 1(1 + x)^-^ - 11 + - (4 - 7)x (1 + x)'-^ + 11 + xlM I . (6) 

Jo X^ L -I J 

Numerically /o(1.65) ~ 25.4, /o(1.8) ~ 33.2, and /o(1.95) ~ 55.6. More specifically the CVT 
predicts the small-scale peculiar velocity dispersion as (Peebles 1976; Suto 1993) 



■"UrVJvr = 1460, 



• /^5.4^ i,^] V'' k,n/.eo. (7) 



^ 1.362 V 33.2 ■ \5Ah-^MpcJ ^l^'^Mpc^ 
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The above simple estimate can be improved more realistically by including the finite 
size of the galaxies, which suppresses the effective gravitational force between pairs, and 
also possibly changes the behavior of the two-point correlation functions on small scales. In 
fact Peebles (1976) already remarked that the above two effects may change the predictions 
of the velocity dispersions significantly. More recently Bartlett & Blanchard (1996) also 
discussed the importance of the local mass distribution around pairs in the CVT. Our 
principal aim here is to describe these effects more extensively using both the perturbation 
analysis and the full numerical integration. For this purpose we take account of the former 
effect by softening the gravitational force according to the Plummer law with the softening 
radius r^. Then equation (|l]) is modified as 



6Gp 



°° dr' , , 
— / az 



(^2 + r2)3/2 



Cp{r\z, \y' - z|). 



We model the latter effect by incorporating a nonzero core radius Tc as follows: 



The two modifications can be still described in equation 



(9) 



simply by replacing /o(7) with 



-?"(?/, 7) 
= 27r(l + y)^ 



(1 + Z/)' 



dt 



(t2 + y^2)3/2 



°° ds 



s-t 



t^dt 



lo {t^ + ^2)3/2 

where y = Vc/r, w = r^/r, and 



ds 



S (|s-t|+|/)T 

1 1 



1 



+ 



1 



(s + yp (t + yp 



;s + ?/)^ ' {t + yp 
Kis,t,y; 7), 



(10) 



jjLdjj, 



1 ( + t2 - 2st/i + yp 
s + t + yf-' -{\s-t\+yf-'^ 



+ 



+ 



+ 





7) ^ 


is 


s^ + e- 


3y' 




2sH^{2- 


-7) 




y{y^ - s2 




) 



2s2t2(l_^) 



s + t + yy-"' -{\s~t\+y) 



2-7 



s + t + yp-^ -{\s-t\+y) 



1-7 



(11) 



The integral I{y,w; 7) is divergent for 7>2if?/ = w = 0. Nonzero values for y and/or 
w make I{y,w;'y) finite for 7 > 2, therefore incorporating the finite size of galaxies in this 
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way widens the applicability of the CVT compared with its simplest form For realistic 
CDM models, the mass two-point correlation functions on small scales may be better 
described by equation (|^) with 7 ;^ 2 and Tc > 0, as will be seen below. 



3. EFFECT OF THE SIZE OF GALAXIES 

3.1. perturbation analysis for 7 < 2 

In the case of 7 < 2, the integral I{y,w;'y) converges even for y = w = 0. Then one 
may treat the effect of nonzero y and w via perturbation theory. Consider the case of w = 
first. Let us set y = (a is a parameter greater than unity and will be determined below), 
and try to find the Taylor expansion with respect to u: 



/(w^0;7) = /(0,0;7)+|^ 

The derivative dl/du is explicitly given as 
1 dl 



u + 0(u^). 



(12) 



u=0 



a-1 



roo roo 

1 + u'^y-^ J dt j ds 



1 



1 
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L *, 
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(S + M")'^ (t + U")T 
1 



ir(s,t,M-;7) 



^(s,t,M";7) 



(s + M")')'+l {t + M")T+1 

K{s,t,u";j + 1). 



1 1 

+ 



{s + u°')-y {t + M")T 

In the limit of m — > 0, only the second term gives a no n- vanishing contribution: 

27 



(13) 



1 dl 

27r du 



roo roo ^ 

— lim a;7M°~^ / dt ds- : 

Jo Ji (t + M" 



7+1 



3(s + u°'p- 



-^t + 0{t' 



(14) 



" 

By setting t = u'^t and then considering the limit m 0, one finds that only the choice of 
= 1/(2 — 7) gives the finite value for equation (|14D which is contributed only by the term 
proportional to t in the square bracket: 

/■°° 2'-)ds rdr 
= -a7 / 

u=0 



1 dl 

27r du 



1 3s7+iyo (r + 1)^+1 3(7-l)(2-7) 



;i<7<2) (15) 
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The above method can be also applied to the perturbative analysis for J(0, w\ 7) and 
/(y, = 7), and a straightforward calculation yields 

l{y,0;^) = /(0,0;7)- 3^^_^^(^_^^ y'-^ + 0{y'-'''), (16) 

J(0,^;7) = /(0,0;7)-^5(^,^)^'"^ + O(«;^-2^), (17) 

/ — 7 z z 



I(y,y;-f) = J(0,0;7)-— -/ -— — ^— f rfr + 0(i/^-^^), (18) 

vy'y- V > , 3(2- 7) io (r2 + l)5/2(r + 1)7+1 vy v ; 

where B{a,l3) is the beta function. The above expressions clearly show the importance of 
the finite size effect for the observational value 7 = 1.8. Since the lowest-order correction 
term begins with y^"''', even a very small value of y changes the value of J(?/,w;7) 
significantly, which will be discussed in details below. 



3.2. numerical evaluation 

Figure Q plots the results of the numerical evaluation of equation (|ToD for 7 = 1.8 in 
comparison with the perturbative expansions ([T6| ) to (18). In practice, here and in what 



follows we set w = 10^^° and y = 10^^°, respectively, when we numerically compute the 
values of /(?/,0;7) and /(0,w;7). The perturbation result is in excellent agreement with 
the numerical evaluation up to y = 0.001 for I{y, 0; 1.8) and I{y, y; 1.8), and up to = 0.1 
for l{0,w; 1.8). Also the perturbation correction turns out to provide a reasonably good 
approximation even up to ?/ ~ w ~ 1. Thus for astrophysically relevant values of and 
Tc, the perturbation formulae (eqs.|]TB[ to |T^) are very useful if 1 < 7 < 2. It should be 
noted that even very small values of y and/or w significantly decrease I{y,w,'~f) relative to 
the idealistic point-mass approximation 1(0,0,7) = -^0(7); for 7 = 1.8, /(y, w, 7)/J(0, 0, 7) 
is about 0.8 even if y or w is ~ 10~^. On the other hand, the dependence on y and w is 
fairly weak as indicated by the perturbation correction (|l^) to ([T8|) . Therefore in a sense 
most of the effect might be interpreted as the inherent offset originated from the point-mass 
approximation itself. 

If 7 > 2, on the other hand, the above perturbative approach does not work and one 
should rely on the fully numerical integration. Figures |^ and ^ show such results; I{y, 0; 7) 
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and I{y,y;'y) for various values of 7. We find the following accurate fitting formulae for 
2 < 7 < 3: 

J(l/,0;7) = A[{y/y,r^ + {y/yay^]-\ (19) 
I{y,yn) = B[iy/y,r + iy/y,f^]-\ (20) 

where 



A 


= 24.99 - 36.20(7- 


2) + 13.76(7- 


2)^ 


(21) 


Va 


= 0.197 + 0.198(7 - 


2) + 0.770(7 - 


2)^ 


(22) 




= 0.0858 + 0.605(7- 


- 2) + 0.274(7 - 


-2)^ 


(23) 


Pa 


= 0.747 + 0.860(7 - 


2) + 0.625(7 - 


2)^ 


(24) 


B 


= 19.23 - 21.65(7 - 


2) + 5.807(7 - 


2)^ 


(25) 


Vb 


= 0.214 + 0.134(7- 


2) + 0.278(7 - 


2)^ 


(26) 


Otb 


= 0.116 + 0.561(7- 


2) + 0.260(7 - 


2)^ 


(27) 




= 0.946 + 1.003(7 - 


2) + 0.469(7 - 


2)2. 


(28) 



At least in the range of 0.001 y ~ 1, the agreement is excellent (Figs. ^ and |^) between 
the numerical results and the fitting formulae. 



4. PREDICTIONS FOR THE SMALL-SCALE VELOCITY DISPERSIONS 

4.1. power-law models 

Let us consider the predictions for the small-scale velocity dispersions on the basis 
of the CVT with the finite size effect. With nonzero and fc, equation (|^) should be 
explicitly rewritten as 



{«L(r)>i?, = 1460 



1.352 ■ ' 1^5.4/i-iMpc^ 



I / \ 
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Due to the r- dependence of I{rc/r,rs/r;'y), {vl2{f))cvT is not described by a power-law in 
general, but is still proportional to ^ilaQpr^/b^. 

Figure I plots (i^?2(^))cvT/^f^o<5p/&g/(ro/5.4/;,-iMpc)^/2 7 = 1.8. As discussed in 
§3, even fairly small values for and Vc can reduce the predictions appreciably relative 
to the case of the point-mass approximation (r^ = Tc = 0). Furthermore, when Vc > 0, 
the softening of the force has negligible effect if < Vc- Since the observed two-point 
correlation functions of galaxies do not show any noticeable departure from the single 
power-law at least for r ^ 0.02/i~^Mpc, the choice of Vg = Tc = 100/i~^kpc plotted here 
might be an extreme example, and already overestimate the finite size effect. 

Similarly the results for various 7 are plotted in Figure ^ for = = 100/i~^kpc 
(dashed curves) and = = lO/i^^kpc (solid curves). The results for larger 7 become 
more sensitive to the values of and r^. Given the fact that the point-mass approximation 
fails for 7 > 2, this is reasonable. It should be also pointed out here that even when 7 > 2 
(f^2('"))cvT increases as r because /(r^/r, r^/r; 7) increases more strongly than oc r'^~^; 
in the case of Vg = Vc = lOO/i^^kpc and r ^ l/i'^Mpc relevant for Figure ^, ('y^2(^))c'vT 
is proportional to r^-^/^+(^b/^ = r°-^^ r°-^\ r°-^^ and r^-^^ ^ ^ 2.2, 2.4, 2.6, and 
2.8, respectively. When r/r^, r/r^ ^ 10, (fi2(^))cvT approaches oc r^"'^/^"'""^'/^ and thus 
significantly levels off (around r ~ l/i^^Mpc for = = 100/i^^kpc and r ~ 0.1/i~^Mpc 
for Tg = Tf. = lO/i'^kpc in Fig.^. In any case, this implies that the scale-dependence of 
(''^12 (^)) is fairly insensitive to the value of 7. 



4.2. COBE normalized CDM models 

The results presented above are fairly general in the sense that they should be valid as 
long as the CVT is correct and equations (§), (H) and (H) are reasonably good approximation 
to the reality. Nevertheless the velocity dispersion (f 12 ('"))cvt is still crucially dependent on 
the parameters including fio, Qpi bg, tq, 7 as well as Vg and Vc which we have discussed in 
details here. Therefore to proceed further one has to assume some sort of specific models. 
For this purpose cold dark matter (CDM) models would provide good working hypotheses. 
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In particular very accurate semi-analytical fitting formulae for ^p(r) have been worked 
out in the literature, which have been successfully tested against high-resolution N-body 
simulations (Hamilton et al. 1991; Jain, Mo & White 1995; Peacock & Dodds 1994, 1996; 
Mo, Jing & Borner 1996). 

CDM models with the primordial scale-invariant (Harrison -Zel'dovich) spectrum have 
the power spectrum P{k) which is asymptotically proportional to for large k. This 
implies that the two-point correlation function in these models should approach constant 
on the corresponding scales; for a purely scale-free primordial spectra of P{k) oc k^, ^p{r) is 
proportional to r~^~" in a linear regime and to in a stable clustering (strongly 

nonlinear) regime. In other words, we naturally expect the small-scale flattening of the 
form d^) even independently of the finite size effect of galaxies. This is another reason why 
the present analysis is very important in realistic cosmological scenarios. 

Given the cosmological parameters Qq, Aq, and h, the normalization constant of the 
fluctuation spectrum can be determined by COBE data (e.g., Sugiyama 1995), and thus 
the (non-linear) mass correlation function ^p(r) can be analytically calculated (Peacock & 
Dodds 1996; Mo et al. 1996). We then determine tq, 7 and Tc by fitting the calculated 
C,p{r) to equation (|]). So determined Vc does not represent the finite size effect of galaxies, 
and if it is below 10/i~^kpc, for instance, we may set Vg = Vc = 10/i~^kpc to include the 
effect. Therefore as long as the hierarchical relation (H) holds for the two- and three-point 
correlation functions of mass, the only remaining parameter is Qp. Although recent 
numerical and semi-analytical studies indicate that Qp is weakly dependent on the scale 
(Suto 1993; Matsubara & Suto 1994; Suto & Matsubara 1994; Jing & Borner 1996; Mo et 
al. 1996), it is about 2 ~ 3 on scales less than l/i^^Mpc. 

We use the following expression for the COBE 2yr data normalization crcoBE(8/i~^Mpc) 
by Nakamura (1996) who performed a fit to the numerical results by Sugiyama (1995): 



(TcoBE(8/i-iMpc) = (3or)V(8r) 




(30) 
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where 



rioh exp 



1 + 



(31) 



no , 

f{x) = (18.1x°-2 + 7.17a;°-^ + 22.1a; + 1.23x^)-\ (32) 

For definiteness we adopt the baryon density parameter Qi, = O.OlS/i"^. Note that the 
COBE 2yr data give a shghtly larger normahzation amphtude (roughly by 10 percent) 
compared with the COBE 4yr data (Sugiyama 1996, in private communication). Our 
conclusions below, however, are insensitive to this issue. 

Figure ^ plots the two-point correlation functions of mass based on the Peacock & 
Dodds (1996) model and our fits to the form (^ for the CDM models in Table 1. The 
parameter fitting is performed over the range 0.01/i~^Mpc < r < 5/i-^Mpc. Although the 
fit does not seem so good especially on small scales, any single power-law form (i.e., Vc = 0) 
yields much poorer fits and the inclusion of the nonzero definitely improve the theoretical 
prediction. 

Figure shows {vl2ir))cvT/^p for the three CDM models in Table 1. The 
model parameters except for Qp are completely determined from theory and the fitting, 
and therefore the result is fairly definite. As mentioned above, recent numerical and 
semi-analytical results suggest Qp = (2 ~ 3) on scales less than lh~^Mpc. With this value 
in mind. Figure |^ implies that Qq = ^ CDM models disagree with the observed values 
(300 ~ 600)km/sec at r = l/i^^Mpc (Davis & Peebles 1983; Mo, Jing & Borner 1993; Fisher 
et al. 1994; Marzke et al. 1995 ) by a factor of (3 ~ 6). Although the conflict with the 

Table 1: Fitted parameters for CDM models 



no 


Ao h 


fTcoBE(8/i-^Mpc) 


7 


ro(/i"^Mpc) 


rc(/i"^Mpc) 


1 


0.5 


1.35 


2.002 


7.67 


0.061 


0.2 


0.7 


0.27 


1.42 


0.79 


0.0 


0.2 


0.8 0.7 


0.78 


1.99 


4.67 


0.087 
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f2o = 1 model has been pointed out frequently in the literature (e.g., Peebles 1976, 1980; 
Davis & Peebles 1983; Suto 1993), our present prediction is the most precise one from a 
theoretical point of view. Finally in Figure H, we plot (fi2(^ = 0.8/i~^Mpc))c'vt/y^Qp as 
a function of both for Aq = and 1 — fio models using the same procedure described 
here. Clearly only the low-density models are compatible with the observations even if one 
allows for a factor of (2 ~ 3) uncertainty of the sample-to-sample variation (Suto 1996; 
Suto & Jing 1996). It is interesting to note from the figure that softening in the force with 
Tg ^ 100/i^^kpc has almost negligible effect on the prediction of (f^g)^^^ three CDM 

models; this is simply because these models have either 7 ^ 2 or Tc ~ lOO/i^^kpc. 



5. DISCUSSION AND CONCLUSIONS 

We have presented a realistic refinement of the cosmic virial theorem by taking into 
account of the finite size effect of galaxies. The CVT in its conventional form which treats 
galaxies as point particles appreciably overestimates small-scale velocity dispersions of 
galaxies. Although this is already remarked by Peebles (1976) and in a slightly different 
context by Bartlett & Blanchard (1996), we presented more detailed analysis on the basis 
of a perturbation theory and systematic numerical calculations. In particular the former 
provides a clear physical insight into the finite size effect. Our perturbation formulae 
(eqs.|]l6[ - |]18|), which are accurate enough [better than 20% in accuracy for J(?/,w;7)] 
for most purposes, can be easily and reliably used to correct for the flattening of mass 
correlation ^p(r) and the force softening for 7 < 2. We have also presented accurate fitting 
formulae for I{y, w; 7) for the case of 7 > 2. Moreover, we have predicted the pairwise 
velocity dispersion for CDM models with the COBE fluctuation normalization. 

As our result indicates that the flattening of mass correlation C,p{r) has significant effect 
on the CVT, one should be very cautious to infer from observational measurements of 
galaxy-galaxy correlation function ^g{r) and galaxy pairwise velocity dispersions. Although 
current observations of galaxy clustering have not revealed any noticeable flattening in 
^g{r) on small scales, non-linear biasing mechanisms may have complicated the relation 
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between C,g{r) and ^p{r). Moreover, the statistical significance of galaxy clustering on very 
small scales r ^ 200/i~^kpc is still rather poor. In fact, intrinsic flattening of correlation 
functions on small scales is expected generally for CDM models with a Harrison -Zel'dovich 
primordial spectrum. We have quantitatively pointed out the importance of taking into 
account such flattening in ^p(r) when one applies the CVT, otherwise the estimated value 
for Qq would be substantially underestimated. 

To avoid the complicated issue of possible non-linear biasing on small scales, one can 
draw constraints on Qq from a comparison of the measured pairwise velocity dispersion 
with our predictions for CDM models of various Qq (See also Mo et al. 1996). As long 
as one relies on the COBE normalization, our CDM predictions are fairly deflnite; for 
realistic CDM C,p{r) with flattening, fortunately, the (unknown) force softening even with 
^ lOO/i^^kpc would not make any difference on the prediction of (^^2)^^^ (^^^ §4-2 and 
Fig.|^). The remaining two parameters which would complicate the comparison between the 
predictions and observations include the normalized amplitude of the three-point correlation 
functions which is unlikely to be determined accurately from a small sample of galaxies, 
and a possible velocity bias of galaxies which may result from the dynamical friction but 
whose quantitative effect is still controversial in literature. 

Although the present analysis provides a theoretical improvement over the conventional 
CVT, the usefulness of the CVT to estimate Qq also relies on the quality of the observational 
data. In fact. Mo et al. (1993) noted that the observed velocity dispersions vary fairly 
signiflcantly from sample to sample, and also are sensitive to the modeling of the data. We 
will discuss this topic in details elsewhere (Suto & Jing 1996). 

We thank Houjun Mo and Gerhard Borner for useful discussions and comments on the 
manuscript, and Takahiro T. Nakamura for providing the fltting formula for the COBE 
2yr normalization. YPJ gratefully acknowledges the receipt of an Alexander-von-Humboldt 
research fellowship. YS is grateful to Gerhard Borner for hospitality at Max- Planck- Inst it ut 
fiir Astrophysik, where the present work was done. This research was supported in part 
by the Grants-in-Aid by the Ministry of Education, Science, Sports and Culture of Japan 
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Fig. 1. — Effect of nonzero y and w on I[y,w]''y) for 7 = 1.8. Symbols indicate the results 
of direct numerical evaluations while thin curves plot those of the perturbation theory. 
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Fig. 2. — /(y, 0;7) for different 7 (r^ = 0). Thick curves indicate for 7 < 2 while symbols 
(direct numerical integration) and thin solid curves (fitting formula) for 7 > 2. 
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Fig. 3. — I{y,y;'j) for different 7 (r^ = r^). Thick; curves indicate 7 < 2 while symbols 
(direct numerical integration) and thin solid curves (fitting formula) for 7 > 2. 
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Fig. 4. — Relative pairwise peculiar velocity dispersions predicted from the CVT (7 = 1.8). 
Different curves correspond to different sets of values for and Tc indicated in the panel. 
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pair separation r (h"' Mpc) 

Fig. 5. — Relative pairwise peculiar velocity dispersions predicted from the CVT for various 
7. Solid and dashed curves correspond to rg — rc — 10/i~^kpc and 100/i"^kpc, respectively. 
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Fig. 6. — Two-point correlation functions of mass for the three CDM models in Table 1 
computed from analytical models (symbols) and the power- law fits with nonzero Tc (curves) ; 
Qo — 1-0 and Aq = 0.0 (open circles and dashed curve), Qq — 0-2 and Aq = 0.0 (filled squares 
and solid curve), and VLq = 0.2 and Aq = 0.8 (open triangles and dot-dashed curve). 
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Fig. 7. — Relative pairwise peculiar velocity dispersions from the CVT as functions of r for 
the three CDM models in Table 1. = 1.0 and Aq — 0.0 (dashed curve), flo — 0.2 and 
Ao = 0.8 (dot-dashed curve), and = 0.2 and Aq = 0.0 (solid curve). Thick curves indicate 
for = 0, while thin curves for = lOO/i^^kpc. 
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Fig. 8. — Relative pairwise peculiar velocity dispersions at r = 0.8/i~^kpc as functions of 
Qq for COBE normalized CDM models {h = 0.7); Aq = models in dashed and Aq = 1 — 
models in solid curves. Thick curves indicate for — 0, while thin curves for = 100/i~^kpc. 



